Change of 9 dependence in 4D SU(iV) gauge theories across the deconfinement 

transition 
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We study the dependence of 4D SU(iV) gauge theories on the topological 6 term at finite tem- 
perature T. We exploit the lattice formulation of the theory, presenting numerical results for the 
expansion of the free energy up to 0(6 ), for N = 3 and N = 6. Our analysis shows that the 8 de- 
pendence drastically changes across the deconfinement transition: the low-T phase is characterized 
by a large- N scaling with 6/N as relevant variable, while in the high-T phase the scaling variable is 
just 6 and the free energy is essentially determined by the instanton-gas approximation. 
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Important physical issues of strong interactions are re- 
lated to the nontrivial dependence of 4D SU(JV) gauge 
theories on the topological parameter 9, which appears 
in the Euclidean Lagrangian as 
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e = jF^(x)F^(x) - i6-^e^ pa F^ v (x)F^(x), (1) 



where 



q(x) 



64tt 2 



-fiupa- 1 jj,. 
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is the topological charge density. The presence of a 
nonzero 9 would break both parity and time reversal 
and the experimental upper bound on it is very small, 
\9\ < 10 -9 Nevertheless, the issue of 9 dependence is 
interesting and relevant to hadron phenomenology, an ex- 
ample being the so-called U^l)^ problem. In the frame- 
work of the large- N expansion [2|t4| , the nontrivial 9 de- 
pendence provides an explanation to the fact that the 
U(1)a symmetry of the classical theory is not realized 
in the hadron spectrum (see, e.g., Rcfs. [5, 6] for recent 
reviews) . 

In this paper we investigate the topological properties 
and the 9 dependence of 4D SU(iV) gauge theories at fi- 
nite temperature T, in particular across the deconfming 
temperature T c . Such properties are known to be rele- 
vant to the thermodynamic behavior of hadronic matter. 
For example, the effective restoration of the U(1)a sym- 
metry in strong interactions at finite T, and in particular 
around the chiral transition, is relevant to the nature of 
the transition itself d, 0, H- 

As we shall better discuss in the following, one expects, 
on general grounds, a crossover between a low-T and a 
high-T regime for 9 dependence, characterized by differ- 
ent large-TV scalings. In particular, the high-T regime 
should be describable by a semiclassical instanton gas 
picture, which instead fails in the low-T regime, where 
6/N turns out to be the relevant large- N scaling variable. 
In 4D SU(iV) gauge theories the deconfming transition is 
first order for N > 3 and gets stronger as N increases. 



It is therefore reasonable to conjecture that the 9 depen- 
dence may sharply change right around T c , where one 
may expect a singular behavior, such as a discontinuity. 

The purpose of our study is to investigate such a 
scenario numerically, presenting results for N — 3 and 
N = 6 to check the N dependence around the deconfine- 
ment transition. The lowest order 0(9 2 ) contribution to 
the free energy, involving just the topological susceptibil- 
ity, has been already investigated around the deconfine- 
ment transition [0|— ll 2[| . However, as we shall discuss in 
detail later on, the study of higher order terms provides 
a more stringent and definite signature for the change of 
the 9 dependence between the two phases. Such a study 
is the main subject of our investigation. 

The finite-T behavior is specified by the free energy 



F(0,T) = -iln j[dA]exp^-J o ' , 



dt I d i xC e 



(3) 

where T is the temperature, V = V/T is the Eu- 
clidean space-time volume, and the gluon field satisfies 
Afj,(l/T,x) = Afj,(0,x). The 9 dependence can be pa- 
rameterized as 

F(9, T) = F(6, T) - T(0, T) = ^ x (T)9 2 s{9, T), (4) 

where x(T) is the topological susceptibility at 9 = 0, 

(Q 2 )e=o 



X = / d 4 x(q(x)q(0)} e=O 



V 



(5) 



and s(#, T) is a dimensionless even function of 9 such that 
s(0,T) = 1. Assuming analyticity at 6 = 0, s(6,T) can 
be expanded as 

b 2 (T)9 2 + b 4 (T)9 4 + ■ ■ ■ , 



,T) = 1 



(6) 



where only even powers of 9 appear. 

At zero temperature, where the free energy coin- 
cides with the ground-state energy, large- TV scaling ar- 
guments 0, applied to the Lagrangian (fT|) indicate 
that the relevant scaling variable [14j is 6 = 6/N, i.e. 



T(9) « N 2 G(9) 



(7) 
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TABLE I: Summary of known T = results for the ratio x/°" 2 
(where a is the string tension) and the first coefficients &2j for 
N — 3,4,6 _17|. The extrapolation of the results for x/°" 2 , 
using the simple Ansatz a + b/N , gives X/(T 2 = 0.022(2) for 
N — > oo. For more complete reviews of results see Refs. 
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x/- 2 


hi 


i>4 


3 


0.028(2) [5] 


-0.026(3) [18-20J 


0.000(1) [18] 


4 


0.0257(10) [211 


-0.013(7) [211 




6 


0.0236(10) [211 


-0.008(4) 


0.001(3) 



as ./V — > oo. Comparing with Eq. Q, this implies the 
large- N behavior 



X = Xoo + 0{N-% b 23 =0(N- 2 J). 



(8) 



We recall that a nonzero value of Xoo is essential to 
provide an explanation to the XJ(1)a problem in the 
large- N limit [!, 0]. The apparent incompatibility of 
Eq. ([7]) with the periodic nature of the topological 9 
term may be solved by a nonanalytic multibranched 9 
dependence of the ground-state energy ]]J Hj], F(9) = 
N 2 min fc H ( id § Lk ) ■ The large- A" scaling © of the 9 de- 
pendence is fully supported by numerical computations 
exploiting the nonperturbative Wilson lattice formula- 
tion of the 4D SU(A r ) gauge theory at T = 0, see, e.g., 
the results reported in Table Q] for N = 3, 4, 6 (see also 
Refs. [5, 6] for recent reviews). This scenario is expected 
to remain stable against sufficiently low temperatures. 

The large-A" scaling (J7|) is not realized by the dilute 
instanton gas approximation. Indeed, at zero tempera- 
ture, instanton calculations fail due to the fact that large 
instantons do not get suppressed. On the other hand, 
temperature acts as an infrared regulator, so that the 
instanton-gas partition function is expected to provide 
an effective approximation of finite-T SXJ(N) gauge the- 
ories at high temperature [15j . high enough to make the 
overlap between instantons negligible. The correspond- 
ing 9 dependence is [HI 



F{6,T) «x(T)(l-co S 0), 
x(T)«T 4 exp[-8^ 2 /5 2 (T)] 



L -N+4 



(9) 
(10) 



using 8ir 2 /g 2 (T) « (ll/3)JVln(T/A) + 0(lnlnT/ln 2 T). 
Therefore, the high-T 9 dependence substantially differs 
from that at T = : the relevant variable for the in- 
stanton gas is just 9, and not 9/N. The instanton-gas 
approximation also shows that xCOj an d therefore the 
instanton density, gets exponentially suppressed in the 
large- A" regime, thus suggesting a rapid decrease of the 
topological activity with increasing N at high T. Since 
the instanton density gets rapidly suppressed in the large- 
N limit, making the probability of instanton overlap neg- 
ligible, the range of validity of the instanton-gas approxi- 
mation is expected to rapidly extend toward smaller and 
smaller temperatures with increasing N. 



The low-T and high-T phases are separated by a first- 
order deconfinement transition which becomes stronger 
with increasing N with T c converging to a finite 
large- A" limit 24(. This suggests the following scenario: 
the crossover between the low-T large-A" scaling 9 de- 
pendence and the high-T instanton-gas 9 dependence, 
respectively given by Eqs. (JTJ) and (|9]), occurs around 
the deconfinement transition, and becomes sharper and 
sharper with increasing N. See, e.g., Refs. Ill. l25l - l27| for 
further discussions of this scenario. 

It is important, at this point, to stress the following: 
even if the instanton-gas prediction, Eq. ©, receives sig- 
nificant corrections as one approaches T c from above, one 
can still conjecture that the phase transition sharply de- 
limits two regimes with a different large- A^ scaling be- 
havior, i.e. that the free energy is a function of 9/N in 
the confined phase and a function of 9 in the deconfincd 
phase. This conjecture is of course more general than the 
instanton gas picture itself. 

The finite-T lattice investigations of the large- A' be- 
havior of x(T) [ll], E3] indicate a nonvanishing large- N 
limit for T < T c , remaining substantially unchanged in 
the low-T phase, from T = up to T c . Across T c a sharp 
change is observed, and x{T) appears largely suppressed 
in the high-T phase T > T c , in qualitative agreement 
with a high-T scenario based on the instanton-gas ap- 
proximation. 

However, to achieve a more stringent check of the ac- 
tual scenario realized in 4D SU(A r ) gauge theories, we 
consider the higher-order terms of the expansion |6]), 
which provide further significant information on the 9 
dependence. Indeed, the expansion coefficients b 2 j are 
expected to scale like N~ 2 ^ if the free energy is a func- 
tion of 9/N and to be N- independent in the instanton-gas 
approximation. The finite-T behavior of such coefficients 
has never been studied numerically until now and is the 
subject of our investigation. In particular, the simple 9 
dependence of Eq. ^ may be observed at much smaller 
T above T c with respect to the asymptotic one-loop be- 
havior (fT0|) of x(T) which is subject to logarithmic cor- 
rections. 

In particular, we aim at clarifying: i) whether a sharp 
change in the large-A^ scaling of b 2 j is observed across 
T c , signalling a change from a 9/N to a 9 dependence of 
the free energy; ii) how rapidly the values of b 2 j above 
T c converge to the instanton gas prediction, Eq. (J9j) , i.e. 



(2j + 2)!' 



3 = 1,2,..., 



(11) 



for the expansion ([S]). These predictions should be com- 
pared to the T = estimates summarized in Table |U It 
has to be stressed that the results (fTTj) depend just on 
the form ([9]) of the free energy and are independent of 
the renormalized coupling costant g(T). 

Due to the nonperturbative nature of the physics of 9 
dependence, quantitative assessments of this issue have 



3 



largely focused on the lattice formulation of the theory, 
using Monte Carlo (MC) simulations. However, the com- 
plex character of the 9 term in the Euclidean QCD La- 
grangian prohibits a direct MC simulation at 6 ^ 0. 
Information on the 9 dependence of physically relevant 
quantities, such as the ground state energy and the spec- 
trum, can be obtained by computing the coefficients of 
the corresponding expansion around 9 = 0. 

We mention that issues related to 9 dependence, par- 
ticularly in the large- TV limit, can also be addressed by 
other approaches, such as AdS/CFT correspondence ap- 
plied to nonsupersymmetric and nonconformal theories, 
see e.g. Refs. 13, 27-2^], and semiclassical approxima- 
tion of compactified gauge theories 

In order to check the change of 9 dependence across 
the deconfinement transition, we numerically compute 
the topological susceptibility and the first few coefficients 
of the expansion ([6]) above T c , for N — 3 and N = 6 to 
check the N dependence. For this purpose we exploit the 
lattice Wilson formulation of SU(iV) gauge theories 



S L = ~ E RcTrlV(x), 



(12) 



where is the standard plaquette operator [32J. The 
coefficients of the expansion around 9 = can be de- 
termined from appropriate zero-momentum correlation 
functions of the topological charge density at 6 = 0. 
These are related to the moments of the 9 = proba- 
bility distribution P{Q) of the topological charge Q and 
parameterize the deviations of P(Q) from a simple Gaus- 
sian behavior. Indeed 13311 . 
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b 4 = 



(Q 2 ) 



b-2 = - 



(Q 4 )-3(Q 2 ) 5 



L t Lf 12(Q 2 ) 
(Q 6 )-15(Q 2 )(Q 4 )+30(Q 2 ) 3 
360<Q 2 ) 



(13) 



(14) 



where \i 1S the the lattice topological susceptibility 
(xi ~ ° 4 X! a is the lattice spacing). The correlation 
functions involving multiple zero-momentum insertions 
of the topological charge density can be defined in a 
nonambiguous, regularization independent way 3J|, and 
therefore the expansion coefficients &2i are well defined 
renormalization-group invariant quantities. This implies 
that they approach their continuum limit with 0(a 2 ) cor- 
rections. 

We evaluate the above quantities in MC simulations for 
several values of the coupling /3 on asymmetric L t x T 3 
lattices (35|. Accurate estimates of b^j require huge 
statistics, because of the large cancellations when eval- 
uated from the expectation values of powers of Q, as in 
Eq. (fTB"]). Therefore, we have to consider a relatively fast 
method to estimate the topological charge Q of a lattice 
configuration. We choose the cooling method, and in 
particular the implementation outlined in Ref. 21|. The 



TABLE II: N = 3 results. We report the value of lattice 
coupling f3, the temporal (L t ) and spatial (L a ) size of the 
lattice, the reduced temperature t = T/T c — 1 [3q| . Xh x/Tc 
and the first two coefficients 62 and bi of the expansion ©. 






Lt, L s 


t 


10 5 Xi 




-12 b 2 


360 64 


6.173 


10, 40 


-0.053(3) 


2.292(7) 


1.84(2) 


0.37(12) 


-4(11) 


6.241 


10, 40 


0.045(3) 


0.645(3) 


0.77(1) 


1.27(7) 


0.7(1.8) 


6.273 


10, 40 


0.095(4) 


0.375(3) 


0.54(1) 


1.15(7) 


1.4(1.4) 


6.305 


10, 40 


0.145(6) 


0.232(2) 


0.40(1) 


1.02(5) 


3.6(7.2) 


6.305 


10, 30 


0.145(6) 


0.233(3) 


0.40(1) 


1.10(7) 


2.9(1.4) 


6.437 


12, 48 


0.147(13) 


0.103(3) 


0.37(2) 


1.07(14) 


-1.1(1.4) 



TABLE III: N — 6 results. We report the same quantities as 
in Table [H] [13 . 



Lt, L s 



10 5 Xi W i x/Tc -12&2 360 bj 



24.797 6, 24 -0.032(10) 17.14(16) 195(8) 0.07(34) -14(18) 

24.912 6, 24 0.045(14) 0.622(13) 9.6(5) 1.15(8) 2.2(1.8) 

24.912 6, 20 0.045(14) 0.631(16) 9.8(6) 1.17(8) 2.2(0.7) 

25.056 6, 24 0.089(8) 0.132(3) 2.41(9) 1.02(4) 1.0(2) 

24.768 5, 20 0.141(7) 0.121(3) 1.28(4) 1.02(2) 1.0(1) 

25.200 6, 24 0.160(8) 0.0316(12) 0.74(3) 1.02(4) 1.1(1) 



stability of the results under cooling is carefully checked; 
we take our data after 15 cooling steps, but differences 
with the results after 10 and 20 cooling steps remain al- 
ways within the errors reported [3(| (for some examples 
see the appendix). Moreover, the stability substantially 
improves with increasing TV, as already noted in the lit- 
erature, also by detailed comparisons with the more rig- 
orous overlap method (which is much more demanding 
numerically), see, e.g., Ref. [f|. 

A summary of the results for N = 3 and N = 6 is pre- 
sented in Tables ITT1 and UTT1 respectively |37j |. The aspect 
ratio L s /Lt in our MC simulations is sufficiently large 
to give rise to infinite- volume results for the observables 
considered within the statistical errors, as shown by the 
comparison of results for different values of L s . For the 
case of both SU(3) and SU(6) we check the continuum 
limit by comparing the results obtained by using two lat- 
tices of different temporal extent at the largest value of 
T, see Tab. lIMnl 

The MC results clearly show a change of regime in the 9 
dependence, from a low-T phase where the susceptibility 
and the coefficients of the 9 expansion vary very little, 
to a high-T phase where the coefficients &2j approach 
the instanton-gas predictions. Fig. [T] shows the data for 
62 ■ In the high-T phase they are definitely not consistent 
with the scaling (|8]), which would imply a factor of four in 
62, in going from TV = 3 to N = 6. On the other hand, in 
the low-T phase 62 does not significantly differ from the 
T = value. This is consistent with the behaviour of \i 
at N = 3, for which we obtain: X i(T = 0.9bT c )/ X i(T = 
0) = 0.98(1) (at f3 = 6.173). A similar behaviour is 
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observed for TV = 6: X i(T = 0.97T c )/ X i(T = 0) = 1.00(2) 
(at $ = 24.797). 

Although our MC results in the high-T phase are ob- 
tained for relatively small reduced temperatures t = 
T/T c — 1, i.e. t < 0.2, the data for 62 show a clear and 
rapid approach to the value &2 = —1/12 of the instanton 
gas model for both TV = 3 and TV = 6, with significant 
deviations visible only for t < 0.1. The high-T values of 
&2 substantially differ from those of the low-T phase, and 
in particular from those at T = reported in Table Q] 
Also the estimates of 64 are consistent with the small 
value &4 = 1/360. 

Our data confirm that \ rapidly decays with increas- 
ing t in both iV = 3, 6 cases. In particular for TV = 6 
we obtain X i(T = 1.09T c )/xz(T = 0) = 0.0136(4) (at 
f3 = 25.056). This suppression is in qualitative agree- 
ment with the one- loop instanton-gas result (TlOl) . but 
larger temperatures are required for a reliable quanti- 
tative comparison, essentially because of the logarithmic 
corrections to Eq. (flQ|). The sharp behavior of the 9 
dependence at the phase transition suggests that T- is 
actually a function of 0/N, as put forward in Ref. [391 ] . 

We have tried to understand the deviations for hi , vis- 
ible at t < 0.1, by taking into account corrections to 
the instanton-gas formula (j9|) through a virial-like ex- 
pansion: the asymptotic formula is corrected by a term 
proportional to the square of the instanton density. For 
example, we may write 



•F(0,T)«x(l 



x 2 ^) + o(x 3 ), 



(15) 



where we use the fact that x(^) is proportional to the in- 
stanton density, and n(0) can be parametrized as n{9) — 
J2k=2 c 2fc sin(6*/2) 2fc . The above formula gives 



12 



1 X 

o c 4 



T 4 



O 



x_ 

T* 



(16) 



If x gets rapidly suppressed in the high-T phase, as sug- 
gested by Eq. (fit)]) and confirmed by the MC results, 
Eq. (|16l) would imply a rapid approach to the asymp- 
totic value of the perfect instanton g shown by 
the data, see Fig. [T] Assuming c 4 weakly dependent 
on TV, Eq. (|16p predicts an exponentially faster conver- 
gence with increasing TV, as also supported by the data. 
Moreover, a hard-core approximation of the instanton 
interactions [22| gives rise to a negative correction, i.e. 
C4 < 0, explaining the approach from below to the per- 
fect instanton-gas value &2 = —1/12. 

In conclusion, our numerical analysis provides strong 
evidence that the 9 dependence of 4D SU(iV) gauge the- 
ory experiences a drastic change across the deconfine- 
ment transition, from a low-T phase characterized by a 
large- TV scaling with 0/N as relevant variable, to a high-T 
phase where this scaling is lost and the free energy is es- 
sentially determined by the instanton-gas approximation, 
which implies an analytic and periodic 9 dependence. 



• n=3, L,=10 

■ N=3, L,=12 

A N=6, L,=6 

v N=6, L,=5 




FIG. 1: (Color online) Finite-T results of 62 for TV = 3 and 
TV = 6, versus the reduced temperature t = T/T c ~ 1. The 
shadowed regions indicate the T = estimates of 62 for TV = 3 
and TV = 6, see Table UJ 



The corresponding crossover around the transition be- 
comes sharper with increasing TV (see Fig. [T]), suggesting 
that the perfect instanton-gas regime sets in just above 
T c at large TV, while %(T) gets drastically suppressed. A 
virial-like expansion suggests that the approach is expo- 
nential in TV; this issue deserves further investigation. 

It is interesting to remark that hints for an early 
onset of an instanton gas regime above the chi- 
ral/deconfinement transition have been provided by re- 
cent MC simulations of full QCD, by looking at the be- 
havior of the relevant susceptibilities [joj . 

Acknowledgements. A large part of the MC simulations 
was performed at the INFN Pisa GRID DATA center, us- 
ing also the cluster CSN4, for a total amount of about 
600 years of CPU time. HP would like to thank the Re- 
search Promotion Foundation of Cyprus for support, and 
INFN, Sezione di Pisa, for the kind hospitality. We thank 
Francesco Bigazzi and Ariel Zhitnitsky for enlightening 
discussions. CB and MD thank the Galileo Galilei In- 
stitute for Theoretical Physics for the hospitality offered 
during the workshop "New Frontiers in Lattice Gauge 
Theories" . 



APPENDIX: STABILITY ANALYSIS OF THE 
TOPOLOGICAL MEASUREMENTS 

In Fig. @ and © the stability of the results presented 
in the paper for 62 and 64 under cooling is shown: n c is 
the number of cooling steps and 62 and 64 are estimated 
by means of (|14|) . using for the determination of topo- 
logical charge Q the prescription of [2l| . Figures refer to 
the case t ~ 0.1, but all the simulations present similar 
behaviour. 
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FIG. 2: Results for SU(3). 
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FIG. 3: Results for SU(6). 
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